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It is typically assumed, without justification, that a weak coupling between a system and a bath 
is a necessary condition for the equivalence of a canonical ensemble and a microcanonical ensemble. 
For instance, in a canonical ensemble, temperature emerges if the system and the bath are uncoupled 
or weakly coupled. We investigate the validity region of this weak coupling approximation, using a 
coupled composite-spin system. Our results show that the spin coupling strength can be as large as 
the level spacing of the system, indicating that the weak coupling approximation has a much wider 
region of validity than usually expected. 
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I. INTRODUCTION 

It is of fundamental importance to revisit statistical 
mechanics, based on principles of quantum mechanics 
(see e.g., [F 8]). Starting from a microcanonical ensem- 
ble, the concept of temperature emerges from a canonical 
ensemble in the limit of weak coupling between a system 
and a bath. This is valid both for classical distributions 
on phase space and for quantum density matrices (see, 
e.g., [l|-y])- A general composite system made of a central 
system and a bath can be described by 






H = H S + H L 



H 



SB, 



(1) 



with Hs the system Hamiltonian, Hb the bath Hamil- 
tonian, and Hsb the system-bath coupling. In order to 
obtain the temperature, the small term Hsb is usually 
neglected. 

In quantum mechanics, the density matrix of a ther- 
mal equilibrium state in a microcanonical ensemble with 
energy between [E, E + S] (5 <C E but much larger than 
the level spacing of the bath) is described by 



" c = -^ 



(2) 



where d is the dimension of the subspace with total en- 
ergy in the interval [E,E + S], and Vc is the projec- 
tion onto that subspace. In the weak coupling limit, 
Hsb — >> 0, the canonical ensemble distribution is derived 
from the microcanonical ensemble by following the stan- 
dard procedure (see, e.g., 0, Q) 



FIG. 1. (Color online) Schematic diagram of a spin system 
(red arrows with circles) coupled to a spin bath (blue arrows). 
The dashed lines denote the coupling between them. 



which states that an arbitrary pure random state in the 
microcanonical subspace is enough to derive the above 
canonical ensemble result (see, e.g., ^H). Several works 
have followed this direction (see, e.g., |9l-[2l|). However, 
the meaning of "weak coupling" has been invoked with- 
out any justification, except in the papers by Dong et 
al. [11] and Reimann [23J]. 

In this paper, we investigate the effect of coupling be- 
tween the system and the bath and justify the validity 
of the weak-coupling approximation in coupled-spin sys- 
tems. These spin systems are very different from the 
harmonic oscillators in Ref. [22[ . This paper is organized 
as follows. We describe the model of a coupled-spin com- 
posite system in Sec. Inland present our numerical results 
in Sec. IIIII A brief discussion and a conclusion are given 
in Sec. ITVl 



s _ 



Tr s (p' 

e -fiH s 

z 



C\ 



(3) 



II. 



COUPLED SPIN SYSTEMS 



with the partition function Z = Tr[exp(— f3Hs)\ and the 
inverse temperature j3 — 1/T. A much stronger state- 
ment, called canonical typicality, has also been made, 



We consider a model for many spin- 1/2 particles (see 
Fig.[T]). Some of these spins are labeled as "system spins" , 
denoted by S m , with m = 1,2,..., M. Others are labeled 



as "bath spins", denoted by Ife, with k = 1,2, ... ,K. 
This composite system is described by the Hamiltonian 
in Eq. ([T]), where 
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(4) 



m=l fc=l 

with A km describing the coupling between the system 
spin S m and the bath spin l k . Here, Q m and w k are 
the Zeeman splittings of the mth system spin and the 
fcth bath spin, respectively. We have set h = 1 and 
ks = 1 for convenience. Note that the system and the 
bath have dimensions of 2 M and 2 K , respectively. The 
system is only experiencing dephasing, because of the 
commutation between the system H$ and the coupling 

Hsb, [Hs,Hsb] — 0. 

The above model [Eqs. ([T]) and (QJ] is exactly solvable. 
Without the system-bath coupling term, the total energy 
E T is the sum of the system energy E s and the bath 
energy E B , i.e., E T = E s + E B . For an eigenstate, the 
eigenvalue of H is 



E, 



E? + Ef 



13 ~ —j ' ^j 
M 



K 






Smz + 2 j U k Ikz , 



m J mz 



(5) 



fe=i 



where s mz and 7fc z are the eigenvalues of S mz and Ikz, 
respectively. The eigenbasis are the usual computational 
basis. Here, i and j label the eigenstates, and range from 
1 to 2 M and from 1 to 2 , respectively. 

With the system-bath coupling term, the eigenvalue of 
an eigenstate becomes 



M 



E, 



E^ 



m J mz 



± 



1 K 
fc=i 



fig+«2 



(6) 



• M 



with B fc = J2 m =i Akm s mz and the corresponding eigen- 
vector is 



(7) 





M K 


IV^±) = (g)k„ 2 )(g)|/ fc ±), 


m=l k=l 


where 


|4+) = cos(0 fc /2) 


/ fa = ^\+sin(^/2) 


|/fc_) = sm(6 k /2) 


4 2 = 5) ~ cos(0 fe /2) 


with cos 9k = 


B k /yJBl+u>l and 
Similarly, i (from 1 


uj k /^/Bl+ul. 



^2 — - - 

Ikz 



l M ) and 
J± (from 1 to 2 K ~ 1 ) label the eigenstates. 
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FIG. 2. (Color online) Canonical relation between ln(Pi) 
and the system energy Ef for various system-bath coupling 
strengths: n = (no coupling, blue dashed line with circles), 
n = 3 (weak coupling, green solid line with crosses), n = 10 
(medium coupling, red solid line with upper triangles), and 
n — 30 (strong coupling, cyan solid line with down trian- 
gles). Here Pi is the probability that the system is in the ith 
eigenstate with energy Ef . For weak and medium couplings 
(n = 0, 3, 10), the curves are almost linear. See also Table [I] 



TABLE I. Inverse temperatures /3 extracted by fitting the 
slope of the curves in Fig. [2] and the corresponding linear 
correlation coefficient r of each fit. 



Coupling strength n 



Inverse temperature /3 
Linear correlation coefficient r 







10 



30 



1.027 0.993 0.875 0.388 
-0.998 -0.998 -0.990 -0.815 



III. 



NUMERICAL RESULTS 



In our calculations, we choose a set of random numbers 
fl m e [0,0.4], w fc G [0,1], A km € [0,0.037]. The number 
of system spins and bath spins are M = 3 and K = 18, 
respectively. For example, here we set the energy shell 
[E,E + 6] by choosing E = -1.722 and 5 = 0.0350. We 
count the total number of composite-system eigenstates 
whose energy is within the given energy shell [E, E + 5] 
and partition this number according to different system 
eigenstates into 2 M numbers. In this way, we obtain the 
probabilities Pi (i = 1, 2, • • • , 8 for M = 3). We adjust 
the coupling strength by multiplying an integer n (e.g., 
n = 0,3,10,30) for all A km . Thus, n = corresponds 
to the case without coupling, n = 3 weak coupling, and 
larger n denotes stronger coupling. 

We plot in Fig. [5] the relation between \n(Pi) and Ef . 



TABLE II. The width of the Gaussian fit of the curves, 00, 
in Fig. [21 Due to the symmetry, only the 4 lowest curves are 
listed. 



System-bath coupling strength 



n = 
n = 3 
n = 10 
n = 30 



Level index 
12 3 4 



1.89 1.89 1.89 1.89 

1.90 1.89 1.89 1.89 
2.06 1.95 1.90 1.94 
3.11 2.37 2.00 2.34 



One would expect a linear relationship if the canonical 
ensemble [Eq. Q] is equivalent to the microcanonical en- 
semble. We clearly see a linear relationship in the cases 
of no coupling (n = 0) and weak coupling (n — 3) as well 
as deviations from a linear relationship in the cases of 
medium coupling (n = 10) and specially for strong cou- 
pling (n — 30). The stronger the coupling between the 
system and the bath spins, the larger the deviation from 
a straight line. By fitting these data with straight lines, 
we can obtain the inverse temperature (3 from each data 
set, i.e., the negative slope of the fitted line. We also ob- 
tain the linear correlation coefficient r, which shows how 
good the linear relationship is, where \r\ = 1 corresponds 
to a straight line. The values of j3 and r are listed in 
Table U for the four cases we present in Fig. [51 Again, we 
find from these values better linear relationships in the 
cases of no coupling and weak coupling. 

From the above numerical results, we observe that the 
canonical ensemble distribution is valid in the cases of 
no coupling and weak coupling, and approximately valid 
for medium coupling. By comparing the level-spacing 
s « 0.1 of the system and the typical coupling strength 
of the medium-coupling-case, Ak m ^$ 0.37, we find that 
the canonical ensemble distribution is a good approxima- 
tion if the coupling strength is less or of the order of the 
level spacing, A^ m < s. In other words, the canonical 
ensemble is equivalent to the microcanonical ensemble if 
the typical coupling terms are less or equal to the level 
spacing of the system. However, the equivalence between 
the two ensembles is broken once Ak m 3> s, as shown in 
the strong coupling case n = 30 in Fig. [2jd) . Compared 
to Tasaki's estimation A]- m <C s [2], our results give a 
much wider range, which significantly extends the previ- 
ous smallness requirement on the system-bath coupling. 

Next we consider how /3 depends on the composite sys- 
tem's parameters. We plot in Fig. [3] the normalized den- 
sity of states (DOS) distribution for each system eigen- 
state. For clarity, we shift each curve up 0.1 in all panels. 
As expected for spin systems, one clearly sees Gaussian 
distributions for the density of states. It is especially in- 
teresting to note that the widths of the Gaussian distribu- 
tions are almost identical in panels EJa)-(c), where one 
finds that the canonical-ensemble distribution is valid. 
We list in Table [IT] the width of the Gaussian fit of the 
curves in Fig. [3J 
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FIG. 3. (Color online) Normalized density of states for 
system-bath coupling strength n — (a), 3 (b), 10 (c), 30 
(d). The red lines show the energy range [E, E + S] for small 
5 < E. Here, E = -1.722. 



A. Weak-coupling limit 

From Fig. EHa)-(c), it is clear that the bath-level dis- 
tribution for each system eigenstate is the same, except 
for a shifted peak position. In the "weak" coupling limit, 
let us assume that the probability of the system to be in 
the ith eigenstate (with the total energy E) is 



Pi(E) = 



1 



ooV 7r 



: exp 



(E - Ei 



(8) 



where <7o is the width of the bath spectrum. Within the 
energy shell [E, E + 5} , the probability becomes 



Pi = 



E+8 



Pi(x)dx ~ Pi(E) 6. 



(9) 



For two system eigenstates, the ratio of probabilities is 



Pi' 



exp 



(E - Eif ~{E- E v 



exp(-PEi) 
exp(-£Ek) 



(10) 

IJ-LUi') 

where (3 R-. -2£/crg, if E > E^E,,. Using Eq. ([10)1 , 
we can estimate the inverse temperature /3 for a given 
E and a bath spectrum width a$. For the zero-coupling 
case, where the bath-spectrum widths are the same for all 
states, we find j3 = 0.964 in the cases shown in Fig.^a)- 
(c). This value of f3 is in good agreement with the first 
three values shown in Table [TJ 



B. Strong-coupling limit 

In the strong-coupling limit, the assumption of <7o be- 
ing the same for all system eigenstates does not hold, 
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FIG. 4. (Color online) Same as Fig. [2] except that the system 
level distribution is now uneven. Note that this implies that 
an uneven level spacing distribution is equally sensitive to the 
system-bath coupling strength. 



TABLE III. Same as Table [I] except that these values are 
extracted from Fig. [4] 



Coupling strength n 



Inverse temperature /3 
Linear correlation coefficient r 







10 



30 



1.325 1.300 1.101 0.734 
-0.998 -0.999 -0.992 -0.915 



as shown in Fig. [U[d). In this limit, it is impossible to 
obtain the (inverse) temperature, since the canonical en- 
semble and the microcanonical ensemble are not equiva- 
lent. From the features of the density of states, we thus 
conclude that the equivalence between a canonical ensem- 
ble and a microcanonical ensemble happens only when 
the system-bath coupling is "weak" , so that the coupling 
does not significantly change the spectrum properties, 
e.g., the spectrum width <7o. 

We would like to remark that we may obtain a negative 
inverse temperature ((3 < 0) if the total energy of the 
composite system is positive (E > 0), according to the 
evaluation of the density of states. 

We now tentatively conjecture that a system with 
evenly-distributed levels is easier to thermalize, and thus 
more sensitive to the coupling strength. With the help of 
our numerics, we check below this conjecture by compar- 
ing the previous almost-even level distribution case with 
the following uneven one. We find that this conjecture is 
invalid. 

We keep all the parameters the same as before except 
Sl m . We then choose another set of fl m , which gener- 
ates unevenly-distributed energy levels for the system 
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FIG. 5. (Color online) Same as Fig. [3] except for a system 
with uneven level distribution. The total energy (the red 
vertical dashed line) is E — -2.214. Note that the DOS is 
not significantly affected by the uneven distribution of the 
system energy levels. 



TABLE IV. The width of the Gaussian fit of the 4 lowest 
curves in Fig. [5] 



System-bath coupling strength 


Level index 
12 3 4 


n = 


1.89 1.89 1.89 1.89 


n = 3 


1.90 1.89 1.89 1.89 


71=10 


2.06 1.95 1.90 1.94 


n = 30 


3.11 2.38 2.00 2.34 



(see Fig. [4j. The results are presented in Figs. [4] and El 
and in Tables ITlTl and ITVl Comparing the two cases, espe- 
cially the linear correlation coefficients r in Table U and 
IIII1 we do not find strong evidence to support the above 
conjecture that systems with evenly-distributed levels are 
more sensitive to the coupling strength between the sys- 
tem and the bath. 



IV. DISCUSSION AND CONCLUSION 

The composite system we calculate is finite, with M = 
3 spins for the system and K = 18 spins for the bath. 
The total number of eigenstates is 2 21 w 2 x 10 6 . But the 
number of eigenstates in the energy shell [E, E + 5] is not 
very large (~ 10 4 ) for the parameters we choose. From 
the numerical point of view, the relative statistical error 
of the data we present in Figs. [5] and 0] is of the order of 
1%. Choosing a larger composite system would reduce 
the statistical error, with the price of more computation 
time. But we believe that our main conclusions would 
still hold qualitatively. 

In summary, we investigate the effect (on the canoni- 
cal ensemble distribution) of the the coupling between a 



spin system and a spin bath. We find that the canonical 
ensemble distribution still holds even if the typical cou- 
pling strength has the same order of magnitude as the 
system's level spacing. This is much larger than what 
was expected |2|. In addition, we observe that the in- 
verse temperature j3 for a composite spin system can be 
expressed as /3 = —2E/ao , with E being the energy in 
the microcanonical ensemble and <7o the width of the den- 
sity of states in the system's eigenstates. 
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